Casimir interaction with an 1/r-dependence 
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We show that in theory it is possible to obtain a Casimir interaction potential that varies with 
distance as 1/r. We achieve this by invoking hypothetical particles having a harmonic oscillator 
interaction potential. The derivation parallels the derivation of the Casimir-Polder interaction be- 
tween atoms in electromagnetism. 
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I. INTRODUCTION 

Johan Diderik van der Waals found by studying the 
equation of state for real gases that there is an attrac- 
tive force between the gas atoms. The force is still there 
even if the atoms have closed electron shells. He was 
awarded the Nobel Price in 1910 for this and related 
work. The van der Waals (vdW) force was found on em- 
pirical grounds and it was not until 1930 that London 1 ^ 
gave a realistic explanation for this force. This force is 
of electromagnetic origin and is caused by fluctuations 
in the electron density within the atoms. It is called a 
dispersion force. This force is long range and has the dis- 



tance dependence, 



-7. 



the interaction potential varies 



as r~ 6 . Later Casimir and Polder^ realized that if the 
distance is large enough the finite speed of light makes 
the force drop off faster, as r~ 8 ; the interaction potential 
varies as r~ 7 . Both in the van der Waals range and in 
the Casimir limit the force is due to fluctuating dipoles. 
For smaller separations also higher order multipole con- 
tributions appear and the distance dependence becomes 
more complex. The interaction potential can be obtained 
as the sum of the zero-point energy of all the electromag- 
netic normal modes^ of the system. One can divide the 
modes into two categories: one that originates in the 
atoms; one that originates in the surrounding vacuum. 
The first (second) type dominates in the vdW (Casimir) 
region. The Casimir-Polder potential is obtained 4 as 
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where ay (u>) , j = 1,2 is the polarizability of atom j. 
This expression is valid in both the vdW and Casimir 
ranges. In the vdW region we have 
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where we in the last step have used the so-called Lon- 
don approximation 1 ^ for the atomic polarizabilities, 
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In the Casimir range 
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Similar forces act between macroscopic objects. Also 
here they are called vdW and Casimir forces. The 
forces between spherical objects behave just as between 
two atoms and the polarizabilities are now the ones for 
spheres. Finite objects, whatever the shape behave as 
spheres if the separation is large enough. If the separa- 
tion is not large enough the distance dependence of the 
force depends on the geometrical shapes of the objects 
and the material the objects are made of. It depends 
on the dispersion of the modes in the objects and in the 
medium surrounding the objects. 

In Casimir's famous work£ from 1948 he studied an 
idealized system of two parallel perfectly reflecting metal 
plates. In this situation there are no metal bulk modes or 
metal surface modes involved at all; only vacuum modes 
contribute to the interaction and force. This means that 
there are no vdW range. The Casimir result is valid at 
all separations, d, between the plates. He found: 



V c (r) = - 
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for the interaction potential and 



F c (r) 
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for the force. In vacuum there are only transverse modes 
with dispersion uj = cq. If the vacuum is replaced by 
a purely dielectric medium there are still only transverse 
modes but their dispersion is different, uj — cq, where c is 
the speed of light in the medium. The only effect on the 
interaction potential and force is that the speed of light in 
vacuum is replaced by the speed of light in the medium. 
However, the purely dielectric medium is an idealization. 
In any real medium there are both transverse and longi- 
tudinal modes and the force becomes more complex. If 
the idealized perfectly reflecting metal plates in Casimir's 
thought experiment is replaced by plates of real metals 
there is a vdW range for smaller distances, smaller than 
approximately c/u> p i, where u p i is the plasma frequency 
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of the metals. In this range the interaction potential 
(force) varies as dT 2 (dT z ). 

A related problem where we clearly see the effect of the 
dispersion of the normal modes of the objects is two par- 
allel 2D (two-dimensional) metal sheets. For a 2D metal 
sheet there are no bulk modes or surface modes. The 
modes are 2D plasmons. These have a dispersion oj ~ y/q, 
This leads to a fractional separation dependence in the 
vdW range^i^; the interaction potential (force) varies as 
d~ 5//2 (<i -7 / 2 ). In the Casimir limit the interaction is the 
same as for two perfectly reflecting plates. If the disper- 
sion of the normal modes of the sheets were different the 
forces would have a different separation dependence in 
the vdW range. 

To summarize, there are different ways to influence the 
separation dependence of the dispersion forces between 
objects. It depends on the shape of the objects, the ma- 
terial the objects are made from and the properties of the 
surrounding medium. In this work we are more drastic 
and study a system where the fundamental interaction 
is not electromagnetic. We want to see if it is possible 
to obtain an interaction potential that varies with dis- 
tance as 1/r, like the fundamental interaction potentials 
in electromagnetism and gravity. In Sec. [II] we intro- 
duce the basic assumptions, the formalism and notation. 
Sec. Hill is devoted to the derivation of the fields from a 
dipole. These results are used as the starting point for 
the derivation of the dispersion potential and force in Sec. 
llVl Finally, in Sec. [V] follows summary and discussions. 



II. THEORETICAL BASIS 



vidually; they are always in pairs. It takes an infinite 
energy to separate them completely. If we assume that 
they are each other's anti-particles the pair will turn into 
two pairs when one tries to separate them past a criti- 
cal distance. This type of interaction leads to composite 
particles that are neutral and have no direct interaction; 
dispersion interactions due to fluctuations are possible. 

To find the effect from fluctuating electric dipoles we 
have to start with the field from a dipole with dipole mo- 
ment p = qd. It is found to be E (r) = —p. Thus the field 
from a dipole is just minus the dipole moment. There are 
no quadrupole or higher order multi-pole contributions as 
opposed to in the ordinary electromagnetic theory. Fur- 
thermore the field has no spatial dependence. This holds 
for any static distribution of charges within the compos- 
ite particle - the field lacks spatial dependence. This 
implicates that the only static field is a dipole field and 
this dipole field is constant throughout all space, inde- 
pendent of the position of the particles. There are no 
other multi-pole fields. This is very encouraging. 

From Einstein's two postulates in special relativity 
follows that there has to be a magnetic field compan- 
ion to the electric field. Furthermore these two fields 
obey the two homogeneous Maxwell's equations, V x E + 

(1/c) <9B J dt = and V B = 0. This means that we may 
introduce scalar and vector potentials, (p and A, respec- 
tively, where B = VxA and E = -V(p - (1/c) dA j dt, 
all in complete analogy with in electromagnetism. 

Now we have all we need for our derivations. We start 
by determining, in the next section, the fields from a 
time-dependent electric dipole. 



The dispersion forces between atoms have an electro- 
magnetic origin where the fundamental interaction is the 
Coulomb interaction between charged particles, electrons 
and protons. We will in this work introduce hypothetical 
particles with a different fundamental interaction poten- 
tial. We make the basic assumptions that this interaction 
travels in vacuum with the speed of light and that Ein- 
stein's two postulates in special relativity holds also for 
this interaction. We will then parallel the derivation of 
the dispersion forces in electromagnetism. 

We have reasons to believe that any fundamental in- 
teraction potential should form closed classical particle 
orbits. There are only two central force fields where all 
bound orbits are closed: one has the interaction poten- 



tial V 



and the other has V 



type is known as a harmonic oscillator potential and is 
our obvious choice of candidate. 

Our notation throughout this work is in analogy with 
the electromagnetic case. We put a tilde above the quan- 
tities to distinguish them from the electromagnetic coun- 
terparts. We assume that the particles have charge ±q~ 
and postulate that the electric field from a charge q at 
the origin is E = qr. This gives the scalar potential the 
form if (r) = —qr 2 /2. This potential is without bounds. 
This means that the particles can never be found indi- 



III. FIELDS FROM A TIME-DEPENDENT 
DIPOLE 

In this section we derive the fields from a time- 
dependent dipole. We perform the derivation along the 
lines used by Heald and Marion^ on a Hertzian dipole at 
the origin. We start with the retarded potentials, 



(p (r, t) = -| / d 3 r'p (r', t - R/c) R 2 ; 
A (r, t) = i / dV J (r',t- R/c) R 2 ; 



R = r — r' 



(6) 



where p and J are the charge and current densities, re- 
spectively. From the relations 



B = V x A 

E = -Vip-±dk/dt 



(7) 



follow 



E(r,t)=f([p\RR-±[p l ]R 2 R- 
v v 

= Ei (r,t)+Ea (r,t)+E a (r,< 
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B(r,t) = /(- J xR§- 
v L J 
= Bi(r,i)+B 2 ( r ,t), 



x R- r 
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where a dot means the time derivative and square brack- 
ets that the function within the brackets is determined 
at retarded times, t — r/c. We have chosen a spherical 
coordinate system with its z-axis along p. We have used 
the relations 

V([p]R 2 )=R 2 V[p] + [p]V(R 2 ); 
V[p]= [p]V(t-\r-r'\/c)=-^[p]R; 



V x 

V x 



= R 2 V x 



V(i? 2 ); 



x V(t- |r-r'|/c) 



x R. 



(9) 

Now, we apply these relations to a Hertzian dipole at the 
origin pointing in the ^-direction. We have 



Jdv' — > Idtz — > = p and R — > r. 



This gives 



E 3 = 
Bi = 
B 2 = 



1 /r 

' 2 
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x r; 
x r. 



(10) 
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We study the two remaining fields on the axis of the 
dipole and in the equatorial plane. From these we get 
the fields in a general point. We start with the axis and 
let At = d/2c, 



(%) = [q(f + At) (z-d/2) 

V / ax 

-q(t' -At) (z + d/2)]z 
= [([q] + [q] d/2c) (z-d/2) 
([q] - [$] d/2c) (z + d/2)] 
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The fields in the equatorial plane are 
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(15) 



Now we get 

Ei = (- [P] + £ [p]) cos 6r + [p] sin 09; 
E 2 = (i[p]-Hi) 2 [^) c os^-i(i)K^' 



! ^) 2 [^]cos0f + i(r) 2 ^] S in( 



2 Vc 

and the total fields are 
E 



(16) 



B 



(r,t) = {-{p}+2^)[i>}-^) 2 [p}}cos9f 

+ {H-Hi)[p]+Hi) 2 [^}sin^; (17) 

(r^) = {-(i)[p] + Hi) 2 [^]}^^. 



Now we have all we need to calculate the dispersion forces 
between the composite particles. This we do in next sec- 
tion. 



IV. DISPERSION POTENTIAL AND FORCE 

We define the polarizability, a, for a composite par- 
ticle through: p = aE. The dispersion interaction be- 
tween two particles, 1 and 2, is found from realizing that 
a polarization of particle 1 gives rise to an electric field 
at the position of particle 2. This field polarizes parti- 
cle 2, which results in a field at the position of particle 
1. Closing this loop results in self-sustained fields - nor- 
mal modes. We find these in analogy with the treatment 
of the two-atom system in the book: Surface Modes in 
Physics^. 

We may reformulate the electric field from a time de- 
pendent dipole as 

E (r, t) = (p-f)f [-p (t - r/c) + 2(t)f>(t- r/c) 
-(L) 2 p( t -r/c)_ 
+ \P — (P • ?) f] [—p (t — r/ c) 
+ l(L)p(t-r/c)-l(^) 2 p(t-r/c) , 

(18) 

and its Fourier transform with respect to time becomes 

E (r, oj) = p (u) e iu}r / c {(p ■ f) f [-1 + 2 (-iur/c) 

+ (u;r/c) 2 +[p- (p-f)r][-l ( 19 ) 
+ \ (-iujr/c) + \ (ujr/c f }. 

On tensor form we may write 

E (r, w) = c tuJT l c j<5 -1 + 2 (-iwr/c) + (uir/cf 

+/3 -1 + i (-iur/c) + i (ur/c) 2 }p(w), 

(20) 

where the tensors are: a = r^r^/r 2 ; (3 = S^—r^r^/r 2 = 

I — a. We let all tensors have double tildes to distinguish 
them from our fields. Things become very simple if we 
choose the third principle axis to point along r. Then we 
have 
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= 0. If we now let 
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|q — 1 + 2 (— iujr/c) + (wr /of 
+ /§[-!+! (-twr/c) + | (wr/c) 2 ] } , 



(21) 



we may set up the self-sustained fields between two po- 
larizable particles as 



p 2 (w) = <5 2 (uj) Ei (r 2 , w) = d 2 (w) 7P1 (w) ; 
Pi (uj) = &! (u>) E 2 (ri , uj) = &i (w) 7p 2 (uj) , 



(22) 



where Ej (r j , w) is the electric field at the position of 
particle j caused by the polarized particle i. Eliminating 
P2 gives 

Pi (uj) = &i (uj) J&2 (uj) 7Pi (w) = &i (uj) a 2 (uj) 7 2 pi (co) , 

(23) 

where in the last step we have assumed that the polariz- 
abilities are isotropic. Now, we have 

Ipi (a;) = (i - at (uj) a 2 (uj) 7 2 ) pi (uj) = 0. (24) 

This system of equations has non-trivial solutions, the 



normal modes, if 



= 0. Now, 
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and 



A(lu) = {1 - <5i (uj) a 2 (uj) e l2wr / c [-1 + 2 (-iuir/c) 
+ (uir/cf] 2 1 x { 1 - a x (uj) a 2 (uj) e i2uir / c 



-1 + \ (-iur/c) + | (ujr/cf 

" ' (26) 

We will need the results on the imaginary frequency 
axis, i.e., 

A(iuj) = (iuj)a 2 (iuj)e- 2 ^ r / c 

4 + 2(|w|r/c)- (ujr/c) 2 
X {l - ai (iuj) a 2 (iui) e - 2 ^ r / c 

-i + HH^-sK/^ ) 

(27) 

Now, the zero-point energy for the system of two inter- 
acting particles, r apart, relative the energy when they 
are at infinite separation is 



(28) 



where the integration is performed in the complex fre- 
quency plane around a contour enclosing the positive part 
of the real frequency axis. This equation is obtained from 
the generalized argument principle^. With a standard 
procedure^ one arrives at an integral along the imaginary 

frequency axis, V (r) = (h/2n) J °° duj In \A(iu) , and the 

resulting interaction potential is 

oo 

V (r) = A / dw In { [l - a x (iui) a 2 (iuj) r 2ur / c 

-14-2 (ujr/c) - (uir/c) 
X [1- ai (iuj)a 2 (iuj)e~ 2ujr /~ c 
_1 + i (wr / c ) _ i (ur/c)^ 

(29) 

Expanding the logarithm, assuming that the polarizabil- 
ities are very small, gives 

oo 

V (r) = -£ j du^ (iui) a 2 (iuj) e ~ 2ujr / c [6 - 12 (ujr/c) 
o 

+ n(ujr/cf -9(ujr/c) 3 + 3(ujr/cf . 

(30) 

Two separation limits emerge, the van der Waals limit 
for small separations, and the Casimir limit for large. In 
the van der Waals limit we have 



V (r) w — (371/271") / dioai (iuj) a 2 (iuj); r^c/ujQ, 

o 

(31) 

where ujq is some characteristic frequency above which 
the polarizabilities are negligible. Note that the in- 
teraction potential lacks an r-dependence in this range 
and consequently there is no force. If we assume 
that the characteristic frequencies are the same for the 
two composite particles and that the so-called London 

approximation^ 2 -, a (iui) = a (0) / 1 + (uj/ujo) 2 , may be 

used for the polarizabilities the potential in the van der 
Waals limit becomes V (r) = — (3/8) a\ (0) a 2 (0) hujQ. 
In the Casimir limit we have 

V (r) = - (25hc/32Trr) & 1 (0) a 2 (0) ; r > c/u , (32) 
and the force is 

F (r) = — (25ftc/327rr 2 ) &x (0) a 2 (0) . (33) 

Here, the potential has the 1/r-dependence that we 
hoped for. The polarizability is unit less and very small. 
A rough estimate of the transition point between the two 
regions is found by equating the Casimir and van der 
Waals potentials, using the London approximation for 
the latter. Doing so one finds that r c « (25c/12itujq). 
Demanding that this is smaller that the closest possible 
distance between the composite particles means that the 
interaction potential may have the 1 / r-dependence at all 
distances. 
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V. SUMMARY AND DISCUSSIONS 

In summary, we have shown that hypothetical parti- 
cles with fundamental interaction potential of harmonic 
oscillator type lead to retarded dispersion potentials be- 
tween composite particles with the same spatial variation 
as the gravitational potential and the Coulomb poten- 
tial. Both these fundamental interactions have serious 
theoretical problems. In the Coulomb case there is the 
problem with the divergent energy in the fields from an 
electron if it is assumed to be a true point particle. The 
present gravity theories have some problems too. One 
very serious problem is that the graviton, the quantum 



of the quantized gravitational field has eluded detection. 
If one would come to the conclusion that gravity is not a 
fundamental interaction after all the present work could 
come handy. 
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